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ABSTRACT 

Aims. In this work we search for signatures of low-dimensional chaos in the temporal behavior of the Kepler-6eld blazar W2R 
1946+42. 

Methods. We use a publicly available, ~ 160000 points long and mostly equally spaced, light curve of W2R 1946+42. We apply the 
correlation integral method to both - real datasets and phase randomized “surrogates”. 

Results. We are not able to confirm the presence of low-dimensional chaos in the light curve. This result, however, still leads to some 
important implications for blazar emission mechanisms, which are discussed. 

Key words. Chaos - BL Lac objects: general, individual: W2R 1946+42 


1. Introduction 

Within the large family of the active galactic nuclei (AGN), 
blazars are the only members whose emission is produced pri¬ 
marily in a relativistic jet via synchrotron and inverse-Compton 
processes. To be identified as a blazar, the jet should be oriented 
at a small angle to the line of sight. The blazar spectral energy 
distribution consists of two peaks - a synchrotron one, gener¬ 
ally peaking in the optical region, and an inverse-Compton one, 
peaking in the X/y-ray region. 

Like most of the other AGN types, blazars are known to be 
highly variable in the optical. In fact, they are the only AGN 
class, for which evidence of signific ant variations clearly ex¬ 
ists even on intra-night time scales dBachev et al.L l2012[ and 
the references therein). Blazar optical variability (i.e. the syn¬ 
chrotron peak variability) on long-term time scales can gener¬ 
ally be described as unpredictable a nd lacking any periodicity 
(see however e.g. iGupta et'aH . l2009ll . In spite of the extensive 
studies throughout the years, no clear reason for such blazar 
variations has been positively identified. Unlike the radio-quiet 
AGN, where accretion disk instabilities are usually invoked to 
account for the variations, for the blazar case similar “variabil¬ 
ity drivers” should be searched for in the processes in a rela¬ 
tivistic jet. They range from purely geometrical ones (changing 
Doppler factor of the emitting blob and/or microlensing from 
intervening foreground objects) to a complex evolution of the 
relativistic particles’ energy distribution (interplay between par¬ 
ticles acceleration and synchrotron/inverse-Compton loses). In 
any case, studying blazar light curves may provide a clue to bet¬ 
ter understand physics of the relativistic jets, provided sufficient 
data are available. 

A way to explore apparently random and non-periodic time 
series is to search for low-dimensional deterministic chaos, 
which may allow to distinguish a purely stochastic process from 


a process driven by a nonlinear (and presumably chaotic) dy¬ 
namical system. While the former case most likely implies the 
presence of many, acting independently factors (e.g. many in¬ 
dependent emitting components), the later one suggests that the 
entire system might be governed by a small (say 3-5) number of 
degrees of freedom. 

Appropriate way to search for low-dimensional chaotic sig¬ 
natures i n a time series is to apply the c orrelation integral (Cl) 
method (iGrassberger & ProcacciaL 1 1 9831 see also §2.2). In the 
past, what concerns AGN, similar analyses have been performed 
on EXOSAT X-ray light curves of sev eral, mostly radio- quiet 
objects, leading to no conclusive results (iLehto et al.Lll993h . due 
perhaps to the insufficient number of data points, ranging from 
~300 to ~4000 for the different sources as well as to the rela¬ 
tively larger photometric uncertainty. 

No conclusive results for the pres ence of low-dimensiona l 
chaos were also obtained for 3C 345 (iProvenzale et all 1 19941) . 
based on the 800 points long optical dataset with linearly 
interpolated missing part. Indications for possible fractal na- 
tur e of the light curves of a few o ther blazars were found 
by iMarchenko & Hagen-ThornI (Il997h : however, again - huge 
missing parts of the light curves were substituted via interpola- 
tio n. Non-linear ( and pe rhaps chaotic) variability has been found 
by iGliozzi et al.l (l2002h in a 10000 points X-ray light curve of 
Akn 564 - a radio-quiet AGN, observed by Ginga satellite. 
Chaotic signature in ~ 20000 points long light curves in the 
X-ray binary GRS 1915+105 (iMisra et all l2004l) was also re¬ 
ported. 

Most of the previous work, as mentioned above, focused on 
Seyfert 1 galaxies, whose observed properties are very different 
than those of BL Lacs (e.g. BL Lacs are very rapidly variable, 
highly polarized, exhibit superluminal motion, emit up to TeV 
energies, while Seyfert 1-s are not) and the physics is also be- 
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W2R 1926+42 - Kepler light curve 



Fig. 1. The long term light curve of W2R 1926+42 from Kepler. 


lieved to be different between them: beamed versus relatively 
isotropic emission. In this paper, we apply the correlation in¬ 
tegral Cl method to search for low-dimensional chaotic signa¬ 
tures in the optical light curve o f W2R 1926+42 — a blazar 
within the Kepler satellite field (lEdelson et al.L l2013l) . By far, 
this object has the largest ever optical dataset (consisting of al¬ 
most 160000 data points in total), due to the unique Kepler capa¬ 
bilities to observe repeatedly the same field. More importantly, 
large segments of equispaced data are available, which is specif¬ 
ically required by the Cl method. In this work, to the best of our 
knowledge, the Cl analysis is applied for the first time on such a 
large equidistant optical dataset of a blazar. Fig. 1 shows the en¬ 
tire lightcurve of W2R 1926+42. Clearly, the object shows rapid 
and significant variability over 5 times in terms of count rates, 
allowing a formal distinction of ’’low” and ’’high” states, which 
will be treated separately later in our analysis. 


2. Search for chaos 

The plan of the paper is the following. In the next section, we de¬ 
scribe the data set under consideration. Subsequently, we discuss 
the method followed to search for chaos in data and results in §3 
and present a discussion in §4. Finally, we end with a summary 
in §5. 


2.1. Data set 


W2R 1 926+42(19:26:31 +42:09:59 , z=0.155) is a BF Fac (an 
“FBF”. IPadovani & Giommil 1995 ) object that h appens to fall 
in the Kenler satellite field (iFdel^n et al.L l2013h . The Kepler 
continuous monitoring covers a time period of about 1.6 years 
(quarters 11 through 17), starting 29.09.2011, with a sampling of 
about 30 minute and a duty cycle of about 90%. The typical pho¬ 
tometric errors are less than 0.002%. From all available datasets 
we identified 18 uniformly sampled segments of length ranging 
from about 1300 to about 16000 points (PDCSAP_FFUX data, 
coveringroughly 4300-8900Arange). Occasionally, where some 
data points (one to rarely upto five) are missing, those are substi¬ 
tuted in terms of linear interpolation, as the Cl method requires 
equidistant sampling. For longer missing parts, the entire dataset 
had to be divided into separate segments. Any remaining parts 
of length less than of 1000 points (if such) are not considered in 
this analysis. 

Although covered by long cadence datasets (“11c”), the short 
cadence data (“sic”, sampling of about a minute) are treated sep¬ 
arately as they provide a much larger time resolution and dataset 






Slope 


Fig. 2. Cl diagrams for random noise (upper panel), Forenz at¬ 
tractor (middle panel) and phase-randomized Forenz attractor 
surrogate (lower panel). For presentation purposes different em¬ 
bedded dimensions, d - 1,2, ....10 (from left to right) are shown 
in different colors. 
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Table 1. Dataset description 


Dataset file name 

Segment 

Start 

End 

N 

kplr006690887-2012004120508.11c .fits 

seg 

1 

1001.23 

1032.27 

1520 


seg 

2 

1033.06 

1063.22 

1477 


seg 

3 

1072.56 

1098.33 

1262 

kplr006690887-2012088054726.llc.fits 

seg 

1 

1126.65 

1153.99 

1339 

kplr006690887-2012179063303.11c .fits 

seg 

1 

1182.76 

1214.98 

1578 


seg 

2 

1215.72 

1244.92 

1439 


seg 

3 

1245.86 

1268.73 

1120 

kplr006690887-2012211050319. sic .fits 

seg 

1 

1274.13 

1289.47 

7510 


seg 

2 

1296.40 

1305.00 

4210 

kplr006690887-2012242122129.sic .fits 

seg 

1 

1306.15 

1336.31 

14760 

kplr006690887-201227712545 3.11c .fits 

seg 

1 

1306.18 

1336.30 

1475 


seg 

2 

1337.20 

1371.32 

1671 

kplr006690887-2012277125453.slc.fits 

seg 

1 

1337.17 

1371.33 

16720 

kplr006690887-2013011073258.llc.fits 

seg 

1 

1373.51 

1404.30 

1508 


seg 

2 

1436.24 

1471.14 

1709 

kplr006690887-2013098041711 .llc.fits 

seg 

1 

1488.69 

1524.96 

1776 


seg 

2 

1526.14 

1557.96 

1558 

kplr006690887-2013131215648.llc.fits 

seg 

1 

1559.25 

1581.58 

1094 


length. For the short cadence datasets where photometric errors 
are larger (up to 0.02%), a median filter is applied over every 
three data points. This procedure helps to largely reduce spuri¬ 
ous effects of significantly deviating single points and to reduce 
almost twice the photometric errors, of course with the price of 
reducing the dataset lengths. The data segments used in this anal¬ 
ysis are given in Table 1. The start and end times are in days 
(BJD-2454833) and N is the number of the points in the seg¬ 
ment. 


2.2. Tests 

The solution of a non-linear dissipative dynamical system is 
known as attractor, if it stays bound within a finite volume of 
the phase space (e.g. Lorenz attractor: lLorenzill963h . An attrac¬ 
tor is strange, if it is of non-integer dimension. The trajectories 
of a strange attractor evolve in a finite volume of the phase space, 
never returning to the same point. The divergence between two 
infinitesimally close trajectories increase exponentially in time 
making the long-term predictions impossible. 

The Cl method is a suitable approach for revealing signa- 
ture of low-dim ensional chaotic behavior in time series data 
(IVio et al.L Il992h . The method relies on the construction of a 
new (empirical) phase space from the available N discrete data 
points. The data set (e.g. the light curve) is separated into seg¬ 
ments of length d. Each segment can be considered as a d- 
dimensional vector (2f;), embedded in the t/-dimensional empiri¬ 
cal phase space. The number of the vector pairs within a distance 
smaller than r is computed as a function of r for different d and 
related to the total number of pairs (n^) for that d. Thus, the cor¬ 
relation integral can be expressed as 


1 ^ 

Q(r) = - y @(r-\Xi-Xj\) 
n„ ^ 




where 0 is the Heaviside's function. Therefore, if the dimension 
of the attractor is D, then 


Cd(r) oc 


f^, d<D 
r^, d>D. 


Hence, increasing the embedded dimension d leads to saturation 
when d > D, which can be used for the estimation of the attractor 
dimension 


Dr = lim 

r^O 


dlogC{r) 

d\ogr 


where Dc is the correlation dimension of the attractor and can be 
a non-integer value. Knowledge of Dc allows the determination 
of the number of differential equations, describing the dynamical 
system, N, which is the first integer value, larger than Dc and 
therefore makes possible drawing conclusions about the physical 
process driving the variability. _ 

The so called Ruelle’s criterion (iRuellel Il9^ relates the 
maximal dimension of an attractor (if present), that can be de¬ 
tected from a time series of N points, with N itself as D,„ax < 
21ogA, which in our case, depending on the segment, means 
Dfim.x “b — 8. 

As pointed out by iLehto et ^ (Il993h . the results from the 
Cl method might be significantly contaminated if a strong peri¬ 
odic component {D - 1) is present in the data. Broadly speak¬ 
ing, false results of chaotic signatures might be an artifact of 
the exact shape of the power-density spectrum (PDS) of the 
dataset. Therefore, one should also perform the Cl test on “sur¬ 
rogate” data, g enerated from the original data PDS after phase - 
randomization (iKaraket al.Ll2010l : ISchreiber & Schmitzlll99^ . 
If the chaotic signature of the original time series is ’’real”, such 
should not be present in the surrogate. 


3. Results 

Before applying the Cl method to real data, we test it with light 
curves of known correlation dimension. Fig. 2 (upper panel) 
shows the results for a 1500 points long random noise light curve 
(infinite dimension). The figure shows the variation of logC/ 
with its slope dXogCIIdXogr, for different embedded dimensions 
id - 1.. 10), each presented by a separate curve (of different color 
in the electronic version). As one sees, each embedded dimen¬ 
sion saturates at its corresponding value, indicating that there 
is no low-dimensional attractor in these data, as expected. Next 
(Fig. 2, middle panel) is for a 1500 points long Forenz attrac¬ 
tor light curve (D^ - 2.06). Here all embedded dimensions of 
d > 2 saturate around 2.06, indeed revealing the presence of 
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a low-dimensional attractor. Finally, in the lower panel of Fig. 
2, we see the results for the same Lorenz light curve, this time 
phase-randomized (surrogate data). As expected, no long dimen¬ 
sional saturation can be seen, meaning that phase randomization 
destroys the attractor (provided such is present, which is the case 
for a Lorenz attractor). 

This simple test suggests that if a low dimensional attrac¬ 
tor is really present in the original dataset, it should not be 
seen in the phase-randomized surrogates, i.e. both logC/ versus 
dlogCIIdlogr plots should differ signihcantly. 

Fig. 3 and 4 show Cl diagrams for different segments of the 
blazar light curve. Although the pattern does not resemble the 
one of the pure random noise, we cannot find any evidence for 
saturation (e.g. the presence of a low-dimensional attractor) in 
any of the segments. As we have tested only for the first 10 em¬ 
bedded dimensions, one can conclude based on this test that the 
correlation dimension (also very similar to so called fractal di¬ 
mension) could indeed be larger than 10. However, the goal usu¬ 
ally is to hnd “low” dimensional chaos with D ~ 3 - 4, which 
may indeed be useful to constrain or develop the theory, as dis¬ 
cussed below. Furthermore, the surrogates show very similar pat¬ 
terns as the real datasets, meaning that any peculiarity should 
be attributed to the particular structures in the light curve PDS. 
Actually, the resu lts are somewhat similar to those reported by 
iLehto et ^(Il993h for their shot-noise model, as well as for most 
of the sources they studied, even though the emission mecha¬ 
nisms in these two (blazar and non-blazar) types of objects are 
very different. 

Finding low-dimensional behavior in a dynamical system (in 
our case - a light curve) can be of huge importance for the the¬ 
oretical models, as it may limit the number of independent vari¬ 
ables (or the number of equations) that governs entirely the sys¬ 
tem. People often invoke the so called single zone model to re¬ 
produce the observed blazar SEDs. In this model the synchrotron 
radiation, which is responsible for the optical emission, is con¬ 
trolled by only a few parameters - electron density, magnetic 
held, Doppler factor, etc. Therefore, if the single zone models 
actually works most of the time, one should expect to hnd a low¬ 
dimensional (say D^3^) behavior in blazar light curves. 

Indeed, the electron energy density evolution due to external 
injection and (synchrotron and synchrotron self-Compton, SSC) 
radia tion loses is descr ibed by the Fokker-Planck (kinetic) equa¬ 
tion (lKardashevlll962ll : 

^ [\ytot\n(,y, 0] + s ( 7 , t) 

where n(y, t) is the relativistic particles energy distribution, 
\ytot\ = \y.sijn\ + Itiicl B^y^ is the total loss and S{y,t) 
is the energy injection (source) term. The synchrotron spec¬ 
trum (atJeastjts_optical_part) is entirely determined by n{y, t) 
(e.g. iTsang & Kirkl l2007l Appendix A). In a case of a single 
monochromatic injection (as frequently assumed) the solution of 
the kinetic equation is ’’trivial” in a sense that it cannot account 
for the intermittent blazar variability (Fig. 1). 

The energy injection (particle acceleration) mechanism is 
generally unknown, however if the injection occurs near the jet 
base it might be considered deterministic in nature, rather than 
completely stochastic. The reason for this assumption is that in 
such a case the accretion flow and/or black hole properties (like 
spin for instance) should control entirely the relativistic particle 
acceleration. Deterministic injection will perhaps lead to deter¬ 
ministic optical variability, evidence for which we did not hnd 
from our results. 


One way to account for the apparently stochastic light 
curves, which would be consistent at a larger extent with our 
hndings, is to invoke a large number of active zones (emitting 
blobs). The overall light curve in this case will be stochastic in 
nature, independently of the exact injection mechanism in each 
zone. 

One such mechanism is based on magnetic reconnections 
that may occur within the jet and accelerate locally electrons 
to relativistic energies, thus leading to intermittent and stochas¬ 
tic in nature variability. The light curve will be perhaps similar 
to the shot noise case - i.e. randomly distributed in time and 
amplitude explosive events (hares) with possibly different char¬ 
acteristic times. The hares can occur completely at random or 
can be a representation of the so called self organized criticality, 
SOC, i.e. when one reconnect ion may lead to an avalanche of 
subsequent reconnections (see lAschwandenL 1201 4 for a recent 
review). The latter possibility is supported by the often reported 
l/f shape of the blazar PDS’s. 

Another way to account for the observed variability is to 
invoke the presence of turbulence in the jet. In these mod- 
els either turbul ent How passes through a standing shock (e.g. 
iMarscheJ. l2014ll leading to instant energy injection in the cor¬ 
responding (shock crossing) plasma volume o r the opposite - a 
shock wave passes th ough turbulent medium dKirk et all Il998l : 
iMarscher et al.L[T992h . The turbulent how in these works is mod¬ 
eled as an ensemble of cells with different (randomly distributed 
within certain limits) properties, which naturally will lead to 
stochastic overall variability. This approach, however, might be 
an oversimplihcation of turbulence. Although no fully developed 
theory of turb ulence exists, tur bulence and chaos seem to be 
deeply related (IProcacciaLll984l) . In fact, low-dimensional frac- 
tal structures have been found in turbulent jets in lab experiments 
dZhao et all l2008h . So, if turbulence is playing a major role in 
shaping blazar jets, one might expect to witness not only frac¬ 
tal spatial structure but perhaps also fractal temporal behavior 
(i.e. low-dimensional chaos in the light curves). On the other 
hand, depending on the exact shock geometry and line of sight 
angle, light c rossing time del ays between emitting cells might 
be expected dMarscheii l2014ll . which may somehow affect the 
overall light curve and eventually - conceal the presence of low¬ 
dimensional chaos. But as we stressed above no fully developed 
turbulence theory is available, so it is difficult to judge whether 
a turbulent blazar jet should manifest chaotic temporal behavior 
or not. 

It is interesting to note that there is no practical difference 
between the Cl diagrams for segments covering the high states 
and the low states of the light curve. One would normally think 
that for the low states, there should be less number (perhaps 
even only one) of active regions, contributing toward the emis¬ 
sion. Should there be many independent regions contributing to 
the emission it would be unlikely to hnd a low dimensional at¬ 
tractor in the light curve. On the other hand, it would be more 
likely to hnd such if a single zone (of a very few) is generating 
the energy. Not hnding a difference in the Cl diagrams means 
there is perhaps no huge difference in the number zones for low 
and high states (i.e. being a signihcant number in both cases). 
Another explanation, of course, is that the temporal evolution of 
a single zone is too stochastic to be described in terms of low¬ 
dimensional chaos. Or may be the entire picture of the blazar 
emission generation is more complex, than we have previously 
thought. 
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4. Summary 

We have applied the Cl method to the longest ever available 
dataset of a light curve of a Kepler blazar - W2R 1946+42 to 
search for signatures of low-dimensional chaotic behavior. Our 
results do not reveal such for any of the segments, covering both 
- high and low states of the blazar. 

Although not entirely conclusive, these results may be sug¬ 
gestive for the mechanism of particle acceleration. Models, 
based on a ’’single engine”, operating perhaps close to the jet 
base seem unlikely. Models, relying on the presence of many in¬ 
dependent active zones, like multiple magnetic reconnections or 
turbulence within the jet are plausible. 

Further tests are however required for a firm conclusion. 
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Slope 


Fig. 3. Cl diagrams for different segments of the blazar light 
curve for different embedding dimensions, same as in Fig. 2. 
The last 3 digits of the dataset (see Table 1) are indicated. The 
upper panel of each box is the real dataset, the lower one - the 
phase-randomized surrogate. None of the diagrams shows clear 
indication for a strange attractor of low dimension, especially 
taking into account that the ’’surrogate” diagrams are very simi¬ 
lar to the ones of the real data. 
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Fig. 4. The same as Fig. 3. 
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